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Abstract 

Purpose:lo assess the association between statins and diverse adverse events in Japanese population. 

Methods: New users of statin who started statin after 6-nnonth period of non-use were identified in 68 hospitals between 
January 2008 and July 2010. In addition to the random sample subcohort, we selected additional subcohort members to 
make the stratified sample subcohort have at least one patient in all subgroups stratified by each combination of statin and 
hospital. By abstraction from medical records, detailed information was obtained for all potential cases and pre-selected 
subcohort members. The event review committee consisting of 3 specialists judged whether possible cases met the 
definition of one of the adverse events of interest, and for adjudicated cases the committee further judged whether statin 
was a certain, probable or possible cause of the occurrence of the event. Adjusted for covariates including age, gender, 
status of "switcher", use of high daily dose and comorbidities at baseline, hazard ratio (HR) was estimated by the Cox 
proportional hazards model with Barlow's weighting method. Data were also analyzed by the method proposed by Breslow 
in 2009. 

Results: A total of 6,877 new users of a statin were identified (median age: 66 years; males: 52%). The hazard ratios of 
increase in serum creatinine for atorvastatin and fluvastatin have wide confidence intervals, but both of the point estimates 
were around 2.5. Estimates of hazard ratios by the method of Barlow (1999) were similar to those by the method of Breslow 
(2009). 

Conclusions: Use of statin was not associated with a significant increased risk for renal, liver and muscle events. However, 
the hazard ratio of increase in serum creatinine tended to be high with atorvastatin and fluvastatin to require further 
studies. 
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Introduction 

Statins are widely used for the treatment of dyslipidemia to 
prevent cardiovascular diseases. In Japan, heart and cerebrovas- 
cular diseases are the second and third causes of mortality [1] and 
the number of patients with hyperlipidemia is estimated as 1.43 
mUlion in 2008 [2] . However, statins have major adverse effects on 
liver and muscles including rhabdomyolysis [3] . Renal toxicity was 
recognized as a concern associated with the use of rosuvastatin 



[4,5] and further addressed recently [6 — 8]. Renal toxicity may be 
augmented with high dose and the maximum daily dose of 
rosuvastatin was set as 20 mg/ day in Japan (half of 40 mg/ day in 
the US and EU) based on a pharmacokinetic study of a small 
number of patients [9]. A company-sponsored post-marketing 
drug use investigation (DUI) [10] was conducted for rosuvastatin 
between 2005 and 2006. As in most DUIs conducted under the 
legislation for re-examination of new drugs, the DUI for 
rosuvastatin did not have a comparator [10]. In this DUI where 
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Table 1. Definition of adverse events. 







Type of event 


Baseline 


Post-dose 


Renal event 


Increase in serum creatinine* 




<4 mg/dL 


>1.5 mg/dL and > 1.5-fold increase from baseline value 




>4 mg/dL 


increase of >0.5 mg/dL 


Hematuria^ 




increase of 2 or more units from baseline on the scale {— ,±, +, ++, +++, mm) 


Proteinuria^ 




at least ++ following increase of 2 or more units from baseline on the scale (—,±, +,++,+++, MM) 


Liver event 


Increase in AST or ALT* 




Normal 


>3 ULN as a result of increase 




Abnormal 


>3 ULN and >2-fold increase when baseline was above the ULN 


Muscle event 


Increase in CK^ 


Normal 


>10 ULN 


Rhabdomyolysis^ 




>10 ULN of CK and clinical course/symptoms 



Abbreviations: ULN, upper limit of normal; AST, aspartate aminotransferase; ALT, alanine aminotransferase; CK, creatinine phosphokinase. 

•Modified from Bellomo et al. 2004 [30]. 

^Criteria used in Klepper MJ and Covert B. 2010 [31]. 

"'"Criteria used in Guidance of FDA. 2009 [32]. Increase of one or bothi of AST and ALT was counted as one event. 
^Criteria used in Pasternak RC et al. 2002 [33]. 
doi:l 0.1 371 /journal.pone.009691 9.t001 



the maximum daily dose was 5 mg/day or less for 98% of patients, 
the increase in serum creatinine by over 30% from the baseline 
was observed in 350 (4.1%) of 8,553 patients during the 12-week 
observation period [1 1]. Though the incidence was unexpectedly 
high, there have been no attempts to compare renal and other 
events between statins in Japanese population. 

We examined the association between statins and multiple 
events (increase of creatinine phosphokinase (CK), rhabdomyol- 
ysis, increase of aspartate aminotransferase (AST), alanine 
aminotransferase (ALT), proteinuria, hematuria and increase of 
serum creatinine) using a prospective stratified case-cohort design. 
A case-cohort design [12,13] has both elficiency, an advantage of a 
case-control study and capability of examining" a variety of events, 
an advantage of a cohort study. 

Methods 

Study Population 

Methods complied with the Strengthening the Reporting of 
Observational Studies in Epidemiology (STROBE) statement [14]. 
We sent a letter of invitation to join the study to a total of 2,037 
hospitals (randomly selected from about 4,000 hospitals "with 150 
or more beds in the nation) and 68 were enrolled. We identified 
patients with hyperlipidemia who newly (after 6-months of non- 
use) started a statin in 68 study hospitals in Japan. Currendy, 6 
statins (pravastatin, atorvastatin, fluvastatin, pitavastatin, rosuvas- 
tatin and simvastatin) are available in Japan and it was in 2005 
that rosuvastatin was marketed as the 6th statin in Japan. In the 
study, first, pharmacists in each hospital identified all patients 
(prevalent users) who used a statin at least once during some time 
window using electronic prescription data maintained inside each 
hospital. The time window was normally a 3-month period but 
could be longer or shorter and each hospital was allowed to select 
any time window provided that it was included in a study period 
between January 1, 2008 and July 31, 2010. Second, to identify 
new users, patients were excluded if prescribed the same statin in a 
prior 6-month period. To know whether the patient was 



prescribed the same statin, we examined both of the electronic 
prescription record and electronic/ non-electronic medical record. 
The latter record was used to exclude patients who already started 
the same statin in a different hospital or clinic before the patient 
was referred to the study hospital. If the patient did not use the 
same statin but used a different statin or other lipid-lowering drugs 
during the preceding 6-month period, the patient was included in 
the final cohort as a "s"witcher" from another lipid-lowering drug. 
Finally, pharmacists sent a list of patients anonymized by the study 
ID number as well as the information on gender, age, generic 
name of statin and the date when the patient used the statin first to 
the research office (NPO Drug Safety Research Unit Japan, 
Tokyo, Japan). 

Sampling of Subcohort 

When the study office received a list of new users from each 
study hospital, 5% of patients in the hospital were randomly 
selected as members of the random sample subcohort. Further- 
more, when no patient was found to be selected as a subcohort 
member by random sampling in one or more of subgroups 
subdivided by statin in each hospital (defined as a "missing 
stratum"), one additional patient was selected (by random 
sampling) from each missing stratum as an additional subcohort 
member. The final stratified sample subcohort consisted of 
random sample subcohort members plus additional subcohort 
members selected from "missing stratum" in each hospital. This 
way of additional sampling would selectively increase the sampling 
fraction from small strata with relatively small additional cost. We 
thought that this feature would be beneficial particularly when one 
of small groups needed to be examined to know whether the group 
had any distinct characteristics while the amount of data obtained 
from random sample subcohort alone would be felt to be too 
small. 

Two-Stage Data Collection Using Standard Report Forms 

We used two types of standard report forms for collecting 
medical records of cohort members. The data was extracted from 
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electronic or non-electronic medical records. The first simple form 
was to obtain the data on all cohort members and ask whether or 
not the patient had blood test for serum CK, AST, ALT and 
creatinine and urine test for hematuria and proteinuria during the 
3-month follow-up period. We also asked whether any of post-dose 
laboratory test results met the criteria of the increase of CK (> 1 0 
upper limit of normal, ULN), AST (>3 ULN), ALT (>3 ULN) or 
serum creatinine (>1.5 mg/dL) as well as hematuria (+) or 
proteinuria {++) during the follow-up period. Potential cases were 
defined as those who met one of these criteria irrespective of the 
baseline laboratory test results. The second form was sent to obtain 
the information on all of potential cases and subcohort members 
pre-selected by the random and additional samplings. The second 
form was to ask diagnoses of and drugs for co-morbidities 
including hypertension, diabetes, heart diseases, liver diseases and 
renal diseases at baseline. We also asked whether or not the patient 
was a "switcher" from another lipid-lowering drug and the date 
when the patient stopped the drug during the observation period 
and date when the patient was lost to follow-up. Specifically for 
potential cases, we collected the data on the clinical course 
(narrative description and detailed laboratory test results before 
and during the follow-up period) by a form for each specific kind 
of adverse event. 

Adjudication of Cases 

We summarized the definition of cases in Table 1 . To identify 
cases likely to be caused by a concomitant condition, the event 
review committee consisting of three specialists in internal 
medicine reviewed all the possible cases while blinded for statin 
used by the patient. A patient was adjudicated as a final case when 
the definition in Table 1 was met and the statin was judged to be a 
certain, probable or potential cause [15]. Patients were excluded 
from the final analysis when a concomitant condition, but not the 
statin, was judged to be a certain or probable cause. 

Sample Size 

To calculate the sample size, we used 0.1% as the lowest 
incidence proportion of AST, ALT and CPK in statin users 
reported in literature [16]. A sample size of the entire cohort 
consisting of two groups of the same size- ^vas estimated to be 
10,436 assuming that the incidence proportion was 0.1% in one 
group and 0.4% or more in another group when a — 0.05, 
power = 0.8 and m = 5 (m was the ratio of subcohort to the 
expected number of cases) and the minimum number of subcohort 
was estimated to be 1 30, a little larger than 1 % of the entire cohort 
[17]. As 6 statins were available in Japan, a required total size of 
the entire cohort was estimated to be 18,000 to 28,000 to make 
comparisons with enough power between the 3 largest subgroups 
when those 3 subgroups consisted of 60 to 90% of all patients. We 
selected 5% of the entire cohort rather than 1% as the random 
sample subcohort because the size of subcohort could be 
substantially smalk-r than the numlx-r of cases for events with 
higher incidence proportion (e.g., proteinuria with 2% of the 
reported incidence proportion) [16]. During the study conduct, it 
was found to be hard to achieve the size of 18,000 or more of the 
entire cohort, but, the attainable sample size (around 7,000) was 
judged to have enough power (around 0.8) to evaluate events with 
higher incidence proportion without altering the fraction of 
subcohort (5%). 

Statistical Analysis 

For the entire cohort, the distribution of age and gender was 
estimated for each of 6 statins and the standardised difference [18] 
was calculated between pravastatin (reference) and other statins for 



the proportion of patients in the entire cohort whose blood and 
urine test results were available during the 3-month follow-up 
period. Standardised differences of less than 0.1 are generally not 
considered meaningful [18]. For subcohort members, the mean 
observation period, co-morbidity at baseline and proportion of 
"switchers" were also estimated. 

The hazard ratio of events was estimated by taking pravastatin 
as a reference. First, the hazard ratio in case-cohort analyses was 
estimated using a Cox regression model with the weighting 
method according to Barlow [19]. The weight was the inverse of 
the sampling fraction of the stratified subcohort or Ni/ni where Ni 
and ni were the size of the entire cohort and that of the stratified 
sample subcohort, respectively, of the i-th statin (i= 1, 2, 6). 
According to Barlow [19], the robust variance was used to 
estimate the 95% confidence interval (CI). To adjust for 
confounding, the following potential confounders were included 
in the Cox regression model: age, gender, status of "switcher", use 
of high daily dose (more than "usual" daily dose recommended in 
the package insert) and concomitant diseases (hypertension, 
diabetes, heart disease, liver disease and renal disease). In addition, 
as an ad hoc analysis, we also estimated the hazard ratios by using 
all the available data in the entire cohort according to Breslow [20] 
assuming that the patient was not lost to foUow-up and observed 
for 91 days for patients who were not a case nor a subcohort 
member. We did not collect the latter data because the Breslow's 
method [20] was published in 2009 after the current study was 
started while the Barlow's method [19] (published in 1999) 
assumed the use of the data for cases and subcohort members only, 
AU the statistical analyses were conducted by SAS (Cary NC, 
USA) except for the analyses by the Breslow's method [20] where 
R software system was used. 

Ethics 

The study protocol was approved by the ethics review 
committee of the Public Health Research Foundation (Tokyo, 
Japan) in April 2007 (No.7C0011). By the ethics committee, die 
waiver of informed consent from individual patients in each study 
hospital was approved. We used anonymized data with serial study 
IDs created by the research ofiice. 

Results 

A total of 1 16,418 prevalent users of any statin were identified in 
the time window set between January 1 , 2008 and July 3 1 , 20 1 0 by 
pharmacists in 68 hospitals. The distribution of 6 statins among 
prevalent users were roughly the same as that of 6 statins marketed 
in Japan in 2008 estimated from the sales amounts, drug prices 
(the price is determined by the government and therefore the same 
throughout the nation) and the "usual" daily dose given in the 
package insert. Of the prevalent users, 6,877 (6%) were identified 
as new users of statins. Table 2 shows data on demographics and 
other characteristics of the entire cohort and subcohort. 
Rosuvastatin had the size (n = 2,050) 14 times larger than 
simvastatin (n=145) in entire cohort. The stratified sample 
subcohort consisted of 8.0% of the entire cohort (551/6,877) 
while the random sample subcohort was 5.1% of the entire cohort 
(348/6,877). The sampling fraction of subcohort was 19.3% with 
simvastatin and 6.9"/o with rosuvastatin indicating that the way of 
additional sampling in the current study selectively increased the 
sampling fraction from small strata. The proportion of users of 
rosuvastatin, atorvastatin and pravastatin was 30%, 26%, and 
2 1 % , respectively of the entire cohort while that of pitavastatin, 
fluvastatin and simvastatin was less than 20%. About 80 (58 — 
90)% of patients had blood test while about 45 (40-53)% of 
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Table 3. Adjudication results of events of interest. 



Classification 


Renal event 






Liver event 


Muscle event 






Increase in 
serum creatinine 


Hematuria 


Proteinuria 


Increase in 
AST or ALT 


Increase in 
CK 


Rhabdomyolysis 


Potential cases* 


351 


661 


387 


343 


63 


11 


Definition of adverse events met^ 


No 


271 


539 


358 


204 


57 


10 


Yes 


80 


122 


29 


139 


6 


1 


Alternative 
cause likely* 


22 


39 


8 


85 


4 


0 


Alternative 
cause (number) 


Renal disease (7), 
Surgical operation (6), 
Other drug (2), 
Urinary tract/bladder 
disease (2), Others (5) 


Indwelling urethral 
catheter* (16), 
Urinary tract/bladder 
disease (1 1), 
Surgical operation 
(4), Others (8) 


Urinary tract/bladder 
disease (3), 
Other drug (2), 
Surgical 
operation (1), 
Others (2) 


Surgical 
operation 
(29), Myocardial 
infarction (21), 
Other drug 
(13), 

Others (22) 


Myocardial 
infarction (3), 
Surgical 
operation (1) 




Adjudicated case 


58 


83 


21 


54 


2 


1 



Abbreviations: AST, aspartate aminotransferase; ALT, alanine aminotransferase; CK, creatinine phosphokinase. 
*Potential cases were those with abnormal post-dose laboratory test results regardless of the baseline level. 
^Definition for baseline and post-dose laboratory test results in Table 1. 

"'"Alternative cause was evaluated for possible cases with adverse events that met definition in Table 1 by the review committee. 
^Insertion of an indwelling urethral catheter. 
doi:1 0.1 371 /journal.pone.009691 9.t003 



patients had urine test at least once during tlie 3-montli follow-up 
period. The standardised difference of the fraction of patients who 
had blood test was —0.13 to 0.07 between pravastatin and the 
other statins except for that between pravastatin and simvastatin 
(0.33). The standardised difference was —0.09 to 0.07 for urine test 
except for that between pravastatin and pitavastatin (—0.20). No 
patient used the daily dose more than the approved maximum 
dose a day. It was 3% or less of patients that used the high daily 
dose (higher than the "usual" daily dose recommended in the 
package insert) except for simvastatin where 1 1 % of patients used 
the high daily dose. The proportion of "switchers" was 16% in 
subcohort or 84% were new users who started a statin after 6- 
month period of non-use of any lipid-lowering drugs. The 
proportion of comorbidity in the subcohort was 32% for diabetes, 
57% for hypertension, 39% for heart disease, and 14% for renal 
disease. 

In Table 3, the results of adjudication of events of interest are 
shown. Of a total of 1,816 possible cases, 1,439 were excluded 
because they did not meet the criteria in Table 1. For 158 of the 
remaining 377, the event was judged to be hkely caused by 
concomitant disease, concurrent drug, surgical operation, insertion 
of an indwelling urethral catheter or other procedure/factor and 
the 219 cases finally adjudicated by the review committee were 
analyzed by the Cox proportional hazards model. 

Table 4 shows the number of adjudicated cases and the 
unadjusted and adjusted hazard ratios for the increase of serum 
creatinine, hematuria, proteinuria and increase in AST/ ALT. 
When pravastatin was used as a reference drug, the adjusted 
hazard ratio of the increase of serum creatinine was between 2.4 
and 2.7 with atorvastatin and fluvastatin but the 95% CI was wide. 
Otherwise no difference was found between pravastatin and other 
statins. The Breslow's method gave estimates similar to those by 
the Barlow's method. As a sensitivity analysis, we analyzed data by 
excluding 16% of "switchers" but the results were essentially the 
same as those in Table 4. 



A 72-year-old man who newly started atoivastatin soon after 
hospitalized for stroke had the increased serum CK from 214 lU/ 
L measured on admission (normal range 56 — 244 lU/L) to 12,954 
lU/L on day 3. The patient used no preceding lipid-lowering 
drugs and had no history of cardiovascular disease. The patient 
also had edaravone and argatroban for the treatment of stroke. 
Atorvastatin and other drugs were discontinued and the serum CK 
level was gradually decreased to 175 lU/L in two months. The 
review committee determined that the patient was likely to have 
developed rhabdomyolysis due to atorvastatin. Another 63-years 
old female patient who had pitavastatin developed the increase in 
serum CK from 364 to 1,724 (normal range 45-163) lU/L in day 
29. Except for these two, no patients had the increase in serum CK 
as defined in Table 1. 

Discussion 

We examined the association between multiple events and a 
variety of statins by the prospective stratified case-cohort design. 
The use of statin had no relationship to the increased risk of our 
targeted events. The adjusted hazard ratio of the increase of serum 
creatinine for atorvastatin and fluvastatin was, however, around 
2.5 when pravastatin was used as a reference drug (Table 4) 
though the 95% CI was wide and inconclusive. 

During the last decade, results of several studies on the renal 
toxicity associated with the use of statin have been published [6 — 
8,21]. In pooled analysis of 30 clinical trials, fluvastatin was 
reported to be safe and effective in chronic renal disease [21]. For 
atorvastatin, a beneficial effect on renal function was reported in 
patients "with diabetes [22] or coronary heart disease [23]. On the 
other hand, in a study using the UK database, the risk for acute 
renal failure was increased "with pravastatin and atorvastatin in 
males and females and fluvastatin in females, compared to no 
statin users [6]. In another study using the Taiwan National 
Health Insurance claims database, the renal risk of atorvastatin 
and rosuvastatin was increased compared to other statins 
(lovastatin, simvastatin, pravastatin and fluvastatin) [7]. In a 
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Table 4. Association between statin and events estimated by the case-cohort analysis using the stratified sample subcohort. 





Event* 


Pravastatin 


Atorvastatin 


Fluvastatin 


Pitavastatin 


Rosuvastatin 


Simvastatin 


Increase in serum creatinine 


Number 


11 


20 




7 




7 




13 




0 




Unadjusted HR (95% CI) 


reference 


1.60 (0.72- 


-3.55) 


2.28 (0.79- 


-6.59) 


1.05 (0.38- 


-2.90) 


0.93 (0.39- 


-2.19) 






Adjusted!" hr (95o/„ cI) 


Barlow's method 


reference 


2.74 (0.95- 


-7.89) 


2.70 (0.66- 


-10.98) 


0.72 (0.21 - 


-2.43) 


0.99 (0.32- 


-3.08) 






Breslow's method* 


reference 


2.38 (0.92- 


-6.17) 


2.68 (0.73- 


-9.79) 


0.71 (0.22- 


-2.29) 


0.84 (0.27- 


-2.58) 






Hematuria 


Number 


18 


12 




4 




20 




28 




1 




Unadjusted HR (95% CI) 


reference 


0.57 (0.26- 


- 1 .27) 


0.79 (0.24- 


-2.60) 


1.85 (0.88- 


-3.87) 


1.22 (0.63- 


-2.39) 


0.62 (0.07- 


-5.43) 


Adjusted' HR (95% CI) 


Barlow's method 


reference 


0.59 (0.25- 


-1.42) 


0.71 (0.20- 


-2.49) 


1.16 (0.53- 


-2.52) 


1.06 (0.51- 


-2.20) 


0.51 (0.05- 


-5.04) 


Breslow's method* 


reference 


0.57 (0.25- 


-1.32) 


0.69 (0.21 - 


-2.29) 


1.13 (0.55- 


-2.32) 


1.02 (0.51- 


-2.04) 


0.48 (0.05- 


-4.40) 


Proteinuria 


Number 


6 


2 




1 




5 




7 




0 




Unadjusted HR (95% CI) 


reference 


0.34 (0.06- 


-1.81) 


0.73 (0.08- 


-6.58) 


1.96 (0.57- 


-6.78) 


1.07 (0.33- 


-3.52) 






Adjusted' HR (95% CI) 


Barlow's method 


reference 


0.35 (0.06- 


- 1 .95) 


0.76 (0.09- 


-6.83) 


1.43 (0.41- 


-5.03) 


0.85 (0.26- 


-2.79) 






Breslow's method* 


reference 


0.37 (0.07- 


-2.03) 


0.81 (0.09- 


-7.09) 


1.35 (0.40- 


-4.52) 


0.98 (0.30- 


-3.15) 






Increase in AST or ALT 


Number 


8 


15 




4 




8 




19 




0 




Unadjusted HR (95% CI) 


reference 


1.64 (0.66- 


-4.09) 


1.79 (0.49- 


-6.51) 


1.65 (0.58- 


-4.69) 


1.86 (0.77- 


-4.48) 






Adjusted ' HR (95% CI) 


Barlow's method 


reference 


1.69 (0.65- 


-4.38) 


1.94 (0.50- 


-7.55) 


1.16 (0.40- 


-3.39) 


1.96 (0.73- 


-5.23) 






Breslow's method* 


reference 


1.75 (0.69- 


-4.39) 


2.07 (0.54- 


-7.92) 


1.22 (0.43- 


-3.49) 


1.94 (0.76- 


-4.96) 







Abbreviations: HR, hazard ratio; CI, confidence interval; AST, aspartate aminotransferase; 
*Definition of event is given in Table 1. 

^Adjusted for age, male, switcher from other lipid-lowering drug, use of high daily dose, 
*Estimates using "standard weights" in ref 20 are shown. 
doi:l 0.1 371 /journal.pone.009691 9.t004 



ALT, alanine aminotransferase. 

hypertension, diabetes, heart disease, liver disease and renal disease. 



recent study on more than 2 million new users of statins using 
several databases in North America and UK, the use of 
rosuvastatin, atorvastatin and simvastatin was associated with the 
increased risk of hospitalization for acute kidney injury compared 
with other statins [8]. Our finding, together with these studies, 
may be regarded as a weak alert for the renal toxicity of 
atorvastatin and fluvastatin. 

Two recent meta-analyses [24,25] suggested that the risk of 
transaminase increase in pravastatin was lower than that in other 
statins. Our findings are not contrary to the results but no firm 
conclusion is possible because 95% CI was wide. 

Muscle-related adverse events (rhabdomyolysis, increase of CK, 
myositis, myalgia and myopathy) [26,27] are of the several safety 
concerns associated with statin use. In our study, only rhabdomy- 
olysis and CK increase were included as muscle-related events of 
interest. Our results suggested that muscle-related adverse events 
were rare as only two cases were adjudicated to have developed 
CK increase including one case of rhabdomyolysis in a total of 
about 7,000 statin new users. 

To address the concern of safety of drugs, the primary data 
collection still has its own role. For instance, claims database 
generally does not include laboratory test results. It is worth 
examining an efficient study design like the case-cohort design to 
broaden options for the studies employable to address a wide 
range of drug-related problems. In particular, Japanese regulatory 



authority has often required drug companies to conduct the DUI, 
a stereotyped cohort study without a comparator group involving 
primary data collection following the registry of the users of a 
newly approved drug. The DUI is still prevailing in Japan without 
thoughtful consideration on whether the study with a comparator 
is a better option. We believe that the case-cohort design can be 
one of the future options to improve the company-sponsored post- 
authorization studies. 

To reduce the bias and enhance the efficiency, our study was 
elaborated in several points. First, to find new users of statins in 
study hospitals, we used a time window set by pharmacists in the 
study hospital. New users were then selected from those patients 
identified in this time window by excluding those who used the 
same statin in a preceding 6-month period. This way to identify 
new users was feasible without any special computational skiU 
using the electronic prescription data. Second, we tried to improve 
the method of subcohort sampling. We expected that the size of 
the subcohort selected by a simple random sampling could be very 
smaU for small subgroups when subdivided by the statin (like a 
subgroup with simvastatin where only 7 patients were selected as 
random sample subcohort members as in Table 2). This might 
create a problem when a small subgroup should be evaluated with 
a special attention. For instance, in a hypothetical situation, the 
incidence of a certain event may be found unexpectedly high in a 
small subgroup (though this did not happen in the current study). 
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We tried a strategy of "one additional subcoliort member from a 
missing stratum in each hospital" and at the end of the study, the 
proportion of a stratified sample subcohort in the entire cohort was 
larger for smaller subgroups. It has been suggested that a case- 
cohort design is subject to a certain t^-pe of the information bias 
[28]. In a case-cohort study, detailed data may be collected from 
subcohort members early in the study while data from cases arising 
outside the subcohort can be collected only after the occurrence of 
an outcome of interest. This may result in the deterioration of the 
information collected from cases outside the subcohort. In our 
study, however, we collected the information from subcohort 
members and all the cases at the same timing (or at the end of 
follow-up) to avoid bias. In our study, this did not present any 
problem as the foUow-up period was relatively short (3-month 
period). Special attention and elaboration of data collection 
methods may be needed in a case-cohort study with a long follow- 
up period. 

In our analysis, estimates by the Breslow's method to use all of 
the available data in the entire cohort [20] did not substantially 
different from those by the Barlow's method using the data of cases 
and subcohort only [19]. This was probably because we obtained 
minimum amount of information except for the cases and 
subcohort members. 

This study has some limitations. New statin use in this study 
might not be perfectly representative of statin use in entire 
country. However, the distribution of prevalent users of statin in 
68 hospitals was almost the same with that in the whole nation 
[29]. Second, pediatric patients (<18 years old) were included in 
this study and the study population may not be homogenous. 
However, the number of pediatric patients was 10 (0.1%) in entire 
cohort (n = 6,877) and 2 (0.4%) in subcohort (n = 551) and the 
effect of the inclusion of pediatric patients on the results would be 
minimal. Third, during the 3-month foUow-up period, the 
proportion of those who had blood test (around 80%) and urine 
test (46%) was less than 100% (Table 2) though the proportion was 
similar between statins. Fourth, some events associated with a drug 
may occur late and the follow-up period longer than 3-month may 
be needed to know the occurrence of such events. However, most 
of renal, liver and muscle events (87% = 854/978) had occurred 
within the 12-week observation period in a DUI conducted in 
Japan [11]. For short-term outcomes like in our study, the 3- 
month follow-up may be therefore aderjuate. Fifth, this study had 
insufficient statistical power to detect the association between use 
of statin and renal, liver and muscle events. The sample size was 

References 

1 . Ministry of Health, Labour and Welfare (MHLW). National vital statisties in 
2008 (in Japanese). A\'ailable; http://wvvvv.mhlw.go.jp/toukei/saikin/hw/ 
jinkou/geppo/nengai08/cll/gaikyou.pclf. Acecssed 2014 February 7. 

2. Ministry of Health, Labour and Welfare (MHLW). Health statistics in 2008 (in 
Japanese). A\'ailable: http:/ /www.mhlw.go.jp/toukei/saikin/hw/hoken/kiso/ 
21.html. Accessed 2014 February 7. 

3. Law M, Rudnicka AR (2006) Statin safety: a systematic review. Am J Cardiol 97: 
52c-60c. 

4. Wolfe SM (2004) Dangers of rosuvastatin identified before and after FDA 
approval. Lancet 363: 2189-2190. 

5. FDA web page on rosuvastatin. Available: http://vmw.fda.gov/ohrms/dockets/ 
dockets/04p0113/04p-0113-pdn0001.pdi?utm_campaign = Google2&utm_ 
source — fdaScarch&utm_mcdium = wcbsitc&utm_term = wolfe crestor 
LDA&titm_content — 1. Accessed 2014 l''ebruary 7. 

6. Hippisley-Cox J, Coupland C (2010) Unintended effects of statins in men and 
women in England and Wales; population based cohort study using the 
QRcscarch database. BMJ 340: c2197. 

7. Chung YH, Lee YC, Chang CH, Lin MS, LinJW, et al. (2013) Statins of high 
versus low cholesterol-lowering efficacy and the development of severe renal 
failure. Pharmacoepidemiol Drug Saf 22: 583-592. 

8. Dormuth CR, Hemmelgarn BR, Paterson JM, James MT, Teare GF, et al. 
(2013) Use of high potency statins and rates of admission for acute kidney injury; 



smaller than initially planned and the 95% CI of the hazard ratio 
was wide for most events though all the cases were carefully 
adjudicated by the event review committee. Lastiy, the inclusion of 

the "switchers" in the study population may have impeded the 
estimation of the correct estimates of the risk as the "switchers" 
may be r(;garded as a non-new user of a drug class of lipid- 
lowering drugs. However, the proportion of the "switchers" was 
relatively small (16%) and the results were essentially the same 
even when excluding the "switchers" in the analysis. 

Conclusion 

In a prospective case-cohort study involving primarv' data 
collection, we found that use of statin was not associated with a 
significant increased risk of renal, liver and muscle events. 
However, a weak alert for renal toxicity of atorvastatin and 
fluvastatin shown in our study require further investigations as our 
study did not have enough power. As a case-cohort study can 
examine the association between a variety of outcomes and 
exposures, the design may be useful in post-marketing studies for 
newly marketed medicines. In addition, the design may b(- useful 
when comparing the subgroups of different size which is often the 
case in the observational study comparing users of drugs in one 
class or drugs with the same indication. 

Supporting information 

Checklist SI STROBE checklist. 

(DOC) 

Aclcnowiedgments 

We thank the Japanese Society lor Pharmaeocpidemiolog)' JSPE) and the 
Japanese Soeiety lor 2"*' eonunitiee on aeademie in Japanese Soeiety of 
Hospital Pharmacists lor cooperation for JSS project. The authors wish to 
express appreciation to the pharmacists of participating hospitals for 
providing data. We also thank Ms. Emiko Shiina for data collection and 
management. 

Autlior Contributions 

Conceived and designed the experiments: NO TS AW TO MK AK HK 
YS KM HY TY SK KK. Performed the experiments: NO TS MK YS SK 
KK. Analyzed the data: NO TS AW TY KK. Wrote the paper: NO TS 
AW TO MK AK HK YS KM HY TY SK KK. 



multicenter, retrospective observational analysis of administrative databases. 

BMJ 346; 1880. 

9. Ministrv of Health, Labour and Welfare (MHLW). Advisory committee on 
applications for approval of drugs (in Japanese). Available; http;/ /www. mhlw. 
gojp/shingi/2004/10/txt/sl018-4.txt. Accessed 2014 February 7. 

10. Tanaka K, Morita Y, Kawabe E, Kubota K (2001) Drug use investigation (DUI) 
and prescription-event monitoring in Japan (J-PEM). Pharmacoepidemiol Drug 
Saf 10: 653-656. 

1 1 . Yoshida S (2007) Safety and effectiveness in drug use investigation on 
Rosuvastatin (in Japanese). Prog Med 27: 1159-1189. 

12. Prentice RL (1986) A case-cohort design for epidemiologic cohort studies and 
disease prevention trials. Biometiika 78; 1-11. 

13. DiPietro NA (2010) Methods in epidemiology'; observational study designs. 
Pharmacotherapy 30: 973-984. 

14. von Elm E, Altman DC, Egger M, Poeoek SJ, Gotzsrhe PC, et al. (2007) The 
Strengthening the Reporting of Observational Studies in Epidemiology 
(STROBE) statement; guidelines for reporting observational studies. PLoS 
Med 4; c296. 

15. Edwards IR, Aronson JK (2000) Adverse drug reactions; definitions, diagnosis, 
and management. Lancet 356; 1255-1259. 

16. Jacobson TA (2006) Statin safety; lessons fi-om new drug applications for 
marketed statins. Am J Cardiol 97: 44C-51C. 



PLOS ONE I www.plosone.org 



7 



May 2014 I Volume 9 | Issue 5 | e96919 



Case-Cohort Study on Multiple Events in Statins 



17. Kubota K, Wakana A (2011) Sample-size formula for case-cohort studies. 
Epidemiol 22: 279. 

18. Mamdani M, Sykora K, Li P, Normand SLT, Strainer DL, et al. (2005) Reader's 
guide to critical appraisal of cohort studies: 2. Assessing potential for 
confounding. BMJ 330: 960-962. 

19. Barlow WE, Ichikawa L, Rosner D, Izumi S (1999) Analysis of case-cohort 
designs. J Clin Epidemiol 52: 1165-1172. 

20. Breslow NE, Lumley T, Ballantyne CM, Chambless LE, Kulich M (2009) Using 
the whole cohort in the analysis of case-cohort data. AmJ Epidemiol 169: 1398— 
1405. 

21. Holdaas H, Wanner C, Ahletshauser G. (jimpelewicz C, IsaaesohnJ (2007) The 
efleet of lluvastatin on cardiac in patients with moderate to severe renal 
insufficiency: a pooled analysis of double-blind randomized trials. IntJ Cardiol 
117:64-74. 

22. Colhoun HM, Bcttcridgc DJ, Durrington PN, Hitman GA, Neil HAW, et al. 
(2009) Effects of atorvastatin on kidney outcomes and cardiovasculeir disease in 
patients with diabetes: an analysis from the collaborative atorvastatin diabetes 
study (CARDS). AmJ Kidney Dis 54: 810-819. 

23. Shepherd J, Kastelein JJP, Bittner V, Deedwania P, Breazna A, et al. (2007) 
Effect of intensive lipid lowering with atorvastatin on renal fiinction in patients 
with coronary heart disease: the treating to new targets (TNT) study. Clin J Am 
Soc Nephrol 2: 1131-1139. 

24. Alberron M, Wu P, Druyts E, Briel M, Mills EJ (2012) Adverse events associated 
with individual statin treatments for cardiovascular disease: an indirect 
compEirison meta-analysis. QJM 105: 145-157. 



25. Naci H, Brugts J, Adcs T (2013) (.^comparative tolerability and harms of 
individual statins: a study-level network meta-analysis of 246,955 participants 
from 135 rgmdomized, controlled trials. Circ Cardiovasc Qual Outcomes 6: 
390-399. 

26. McClure DL, Valuck RJ, Glanz M, Hokanson JE (2007) Systematic review and 
meta-analysis of clinically relevant adverse events from HMG CoA reductase 
inhibitor trials worldwide from 1982 to present. Pharmaeoepidemiol Drug Saf 
16: 132-143. 

27. Shannon JA, John SM, Parihar HS, AUcn SN, Eerrara JJ (2013) A Chnical 
Review of Statin-Associated Myopathy. JoumEd of Pharmacy Technology 29: 

219-230. 

28. Langholz B, Thomas DC (1990) Nested case-control and case-cohort methods of 
sampling from a cohort: a critical comparison. AmJ Epidemiol 131: 169—176. 

29. Anonymous (2011) Yakuji Handbook 2011 (in Japanese). Tokyo: Jiho. 134 p. 

30. Bellomo R, Ronco C, Kellum JA, Mehta RL, Palevsky P, et al. (2004) Acute 
renal failure -definition, outcome measures, animal models, fluid therapy and 
information technology needs: the second international consensus conference of 
the acute dialysis quality intiative (ADQI) group. Crit Care 8: R204-R212. 

31. Klcpper MJ, Cobert B (2010) Drug Safety Data: how to analyze, summarize, 
and interpret to determine risk. Jones&Bardett Learning: MA. 293 p. 

32. FDA. Guidance for industry drug-induced Hver injury: premarketing clinical 
evaluation. Available: http://www.fda.gov/downloads/Drugs/. ../Guidances/ 
UCM174090.pdf Accessed 2014 Eebruary 7. 

33. Pasternak RC, Grundy SM, Smith S, Cleemanjl, Bairey-Merz CN, et al. (2002) 
ACC/AHA/NHLBI Clinical advisory on the use and safety of statins. J Am Coll 
Cardiol 40: 567-571. 



PLOS ONE I www.plosone.org 



8 



May 2014 I Volume 9 | Issue 5 | e96919 



